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Proton-neutron quasi-particle random phase approximation (pn-QRPA) theory has recently being 
used for calculation of stellar weak interaction rates of /p-shell nuclide with success. Neutrino losses 
from proto-neutron stars play a pivotal role to decide if these stars would be crushed into black 
holes or explode as supernovae. The product of abundance and positron capture rates on 55 Co 
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is substantial and as such can play a role in fine tuning of input parameters of simulation codes 
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specially in the presupernova evolution. Recently we introduced our calculation of capture rates 
on 55 Co, in a luxurious model space of 7hoj, employing the pn-QRPA theory with a separable 
interaction. Simulators, however, may require these rates on a fine scale. Here we present for the 
O ' first time an expanded calculation of the neutrino energy loss rates and positron capture rates on 

55 Co on an extensive temperature-density scale. These type of scale is appropriate for interpolation 
purposes and of greater utility for simulation codes. The pn-QRPA calculated neutrino energy loss 
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' rates are enhanced roughly up to two orders of magnitude compared with the large-scale shell model 
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qq , calculations and favor a lower entropy for the core of massive stars. 
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I. INTRODUCTION 
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Supernovae are nature's grandest explosions. They are also responsible for synthesizing most of the elements 
of nature, including those that form our own planet, Earth. Many exotic states of matter, including black holes 
and neutron stars, also owe their existence to supernovae. Since 1934, when Baade and Zwicky [1] suggested that 
supernovae are energized by the collapse of an ordinary star to a neutron star, scientists started debating about the 
physical mechanism responsible for these spectacular explosions (where the luminosity of the star becomes comparable 
to that of an entire galaxy containing around 10 11 stars!). Whereas gravity remains the undisputed source of energy, 
the relative roles of other physical phenomena including, but not limited to, roles of neutrino, continue to be argued. 
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The evolution of a massive star (of masses above 10M Q ) comprises of several stages including hydrogen, helium, 
carbon, neon, oxygen and silicon burning. The time scales, temperatures, densities, fuels, luminosities and neutrino 
losses for these stages can be found in Ref. [2] . During the late phases of evolution of these massive stars an iron core 
develops (of mass around 1.5M Q ). The binding energy per nucleon curve prohibits any further production of energy 
by nuclear fusion, yet the neutrino losses continue unabated, exceeding the Sun's luminosity by a factor of about 
10 15 . Capture rates and photodisintegration processes contribute in the lowering of the degeneracy pressure required 
to counter the enormous self-gravity force of the star. Under such extreme thermodynamic conditions, neutrinos are 
produced in abundance. Eventually the collapse of the iron core begins. The mechanism of core-collapse supernovae 
is strongly believed to depend upon the transfer of energy from the inner core to the outer mantle of the iron core. 
Neutrinos seem to be the mediators of this energy transfer. The shock wave, produced as a result, stalls due to 
photodisintegration and neutrino energy losses. Once again the part played by neutrinos in this scenario is far from 
being completely understood. Whereas Haxton [3] proposed the mechanism of " preheating" by neutrinos as a means of 
assistance for shock revival, Bruenn and Haxton [4] later discouraged the preheating mechanism. They worked on two 
different models, simulating weak and strong shock cases, and found out that in neither case is the energy transferred 
to the matter by neutrino-nucleus absorption significant in terms of preheating the infalling iron-like material. More 
recently, Langanke and collaborators [5] had some observations on the models of Bruenn and Haxton [4] and reported 
much larger preshock heating rates, albeit acting for too short a time to lead to consequences for shock propagation. 
According to authors in Ref. [5] the inelastic neutrino-nucleus scattering modifies the radiated neutrino spectra and 
strongly reduces the high-energy spectral tail of the electron neutrino burst at shock breakout. 

A few milliseconds after the bounce, the proto-neutron star accretes mass at a few tenths of solar mass per second. 
This accretion, if continued even for one second, can change the ultimate fate of the collapsing core resulting into a 
black hole. Neutrinos are the main characters in this play and radiate around 10% of the rest mass converting the 
star to a neutron star (e.g. [2]). Despite the small neutrino-nucleus cross sections, the neutrinos flux generated by 
the cooling of a neutron star can produce a number of nuclear transmutations as it passes the onion-like structured 
envelope surrounding the neutron star. Within <~ 0.1 s of the beginning of the collapse the nonthermal neutrino 
emission is dominated by electron neutrinos owing to the decay and capture of leptons by nuclei and free protons. 
The mean individual neutrino energies are some ~ 10 MeV and constitute roughly 10% of the total available energy 
of around (3 - 5)x 10 53 erg [6]. 

The neutrino energy loss rates are important input parameters in multi-dimensional simulations of the contract- 
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ing proto-neutron star. The reenergizing by charged-current electron neutrino and antineutrino absorption on the 
dissociation-liberated protons and neutrons in the postshock flow remains integral to the supernova paradigm and 
neutrino transport is arguably the single most important component of any supernova model [7]. (For a review of 
supernova neutrino microphysics see also Ref. [8].) Parameter-free multi-dimensional models, with neutrino transport 
included consistently throughout the entire mass, yield conflicting results on the key issue of whether the star actually 
explodes. Reliable and microscopic calculations of neutrino loss rates and capture rates can contribute effectively in 
the final outcome of these simulations on world's fastest supercomputers. 

During core infall electron neutrinos are produced almost entirely by electron captures on free protons and nuclei 
and at sufficiently high temperatures antineutrinos are also produced as a result of positron captures on neutrons. 
Electron capture on protons and positron capture on neutrons also play a crucial role in the evolution of star and 
supernova explosion. During the collapse and accretion phases, they decrease the degenerate pressure in the stellar 
core. The neutrinos produced in these capture processes carry the energy away and result in the lowering of the 
entropy of the core. Positron captures are of great importance in high temperature and low density locations. In 
such conditions, a rather high concentration of positrons can be reached from e~ + e + «-» 7 + 7 equilibrium favoring 
the e~e + pairs. Positron capture on elements lying at the bottom of the valley of nuclear beta-stability (so-called s 
elements) may capture positrons and be transformed into a proton-rich isobars (so-called p elements). The electron 
capture on proton and the positron capture on neutron are considered important ingredients in the modelling of 
Type-II supernovae [9]. 

Fuller, Fowler, and Newman (FFN) [10] performed the first-ever extensive calculation of stellar weak rates including 
the capture rates, neutrino energy loss rates and decay rates for a wide density and temperature domain. They made 
this detailed calculations for 226 nuclei in the mass range 21 < A < 60. They also stressed the importance of the 
Gamow- Teller (GT) giant resonance strength in the capture of the electron and estimated the GT centroids using 
zeroth-order (Ohw ) shell model. Later, Aufderheide et al. [11] extended the FFN work for heavier nuclei with A > 
60. They tabulated the 90 top electron capture nuclei averaged throughout the stellar trajectory for 0.40 < Ye < 0.5 
(see Table. 25 therein). Since then theoretical efforts were concentrated on the microscopic calculations of capture 
rates of iron- regime nuclide. Large-scale shell model (e.g. [12]) and the proton- neutron quasiparticle random phase 
approximation theory (pn-QRPA) (e.g. [13]) were used extensively and with relative success for the microscopic 
calculation of stellar capture rates and neutrino energy losses. Monte Carlo shell-model is an alternative to the 
diagonalization method and allows calculation of nuclear properties as thermal averages (e.g. [14]). However it does 
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not allow for detailed nuclear spectroscopy. 

Nabi and Klapdor [15] calculated weak interaction rates for 709 nuclei with A = 18 to 100 in stellar matter using the 
pn-QRPA theory. These included capture rates, decay rates, neutrino energy loss rates, probabilities of beta-delayed 
particle emissions and energy rate of these particle emissions. Since then these calculations were further refined with 
use of more efficient algorithms, incorporation of latest data from mass compilations and experimental values, and 
fine-tuning of model parameters [16, 17, 18, 19]. 

55 Co is abundant in the presupernova conditions and is thought to contribute effectively in the dynamics of presu- 
pernova evolution. Aufderheide and collaborators [11] placed 55 Co among the list of top ten most important capture 
nuclei during the presupernova evolution. Later Heger et al. [20] also identified 55 Co as the most important nuclide 
for capture purposes for massive stars (25M Q ). Realizing the importance of 55 Co in astrophysical environments, Nabi, 
Rahman and Sajjad [21] reported the calculation of electron and positron capture rates on 55 Co using the pn-QRPA 
theory (see also Ref. [22]). However there was a need to perform a fine calculation of these important capture 
rates on a detailed temperature-density grid suitable for collapse simulation codes (see, for example, Ref. [18, 23]). 
Due to the extreme conditions prevailing in these scenarios, interpolation of calculated rates within large intervals of 
temperature-density points posed some uncertainty in the values of capture rates for collapse simulators. Further, as 
mentioned above, the neutrino energy loss rates needed to be included in these expanded calculations on a detailed 
stellar temperature-density grid to make them more useful in simulation codes. 

In this paper we present for the first time an expanded calculation of (anti)ncutrino energy loss rates and positron 
capture rates on 55 Co at fine intervals of temperature-density intervals. Section II deals with the formalism of our 
calculation. Due to existing physical situation and lack of experimental data the uncertainties present in stellar rate 
calculations are considerable. We discuss the uncertainties of the pn-QRPA model in Section III. In Section IV we 
will be presenting some of our results. Comparisons with earlier calculations are also included in this section. We 
finally will be concluding in Section V and at the end Table V presents our expanded calculation of (anti)neutrino 
energy loss rates and positron capture rates on 55 Co. 

II. THE QUASI-PARTICLE RANDOM PHASE APPROXIMATION WITH A SEPARABLE 

INTERACTION 

The QRPA theory is an efficient way to generate GT strength distributions. These strength distributions constitute 
a primary and non-trivial contribution to the calculation of positron capture and neutrino energy loss rates. Kar 
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ct al. [24] pointed out that the quasiparticle random phase approximation (QRPA) method is quite successful in 
predicting the weak interaction rates of ground states all over the periodic table and also stressed the need to extend 
these methods to non-zero temperature domains relevant to presupernova and supernova conditions. QRPA is also 
the method of choice in dealing heavy nuclei [25]. The QRPA theory considers the residual correlations among the 
nucleons via one particle one hole (lp-lh) excitations in a large model space. Nabi and Klapdor [15] extended the 
QRPA model to configurations more complex than lp-lh. 

We used the pn-QRPA theory to calculate the GT strength functions and the associated capture and neutrino 
energy loss rates for 55 Co. The reliability of the pn-QRPA calculations was discussed in detail by Nabi and Klapdor 
[13]. There the authors compared the measured data (half lives and B(GT) strength) of thousands of nuclide with 
the pn-QRPA calculations and got fairly good comparison. We incorporated experimental data wherever available 
to further strengthen the reliability of our calculated rates. The calculated excitation energies (along with their 
log ft values) were replaced with the experimental ones when they were within 0.5 MeV of each other. Missing 
measured states were inserted and inverse and mirror transitions were also taken into account. We did not replace 
the theoretical levels with the experimental ones beyond the excitation energy for which experimental compilations 
had no definite spin and/or parity assignment (2.98 MeV in case of 55 Co). The pn-QRPA theory was used with a 
separable interaction which granted us the liberty of performing the calculations in a much larger single-particle basis 
than a general interaction. We performed the pn-QRPA calculations using a model space of seven major harmonic 
oscillator shells (7huj). The Hamiltonian for our calculations was of the form 

H QRPA = H sp + ypair + (1) 

here H sp is the single-particle Hamiltonian, V pair is the pairing force, Vq^ is the particle-hole (ph) Gamow- Teller force, 
and Vq^ is the particle-particle (pp) Gamow- Teller force. Single particle energies and wave functions were calculated 
in the Nilsson model, which takes into account nuclear deformations. Pairing was treated in the BCS approximation. 
The proton-neutron residual interactions occurred as particle-hole and particle-particle interaction. The interactions 
were given separable form and were characterized by two interaction constants x an d respectively. The selections 
of these two constants were done in an optimal fashion. For details of the fine tuning of the Gamow- Teller strength 
parameters, we refer to Ref. [26, 27]. In this work, we took the values of \ = 0.2MeV and k = 0.07 MeV. Other 
parameters required for the calculation of weak rates are the Nilsson potential parameters, the deformation, the 
pairing gaps, and the Q- value of the reaction. Nilsson-potential parameters were taken from Ref. [28] and the Nilsson 
oscillator constant was chosen as ftw = AlA^ 1 ^ 3 (MeV) (the same for protons and neutrons). The calculated half-lives 



depend only weakly on the values of the pairing gaps [29]. Thus, the traditional choice of A p = A„ = 12/\/A(MeV) 
was applied in the present work. The deformation parameter for 55 Co, 5, was taken to be 0.06, according to Mollcr 
and Nix [30]. (See also the discussion on choice of deformation parameter in Ref. [19].) Q- values were taken from the 
recent mass compilation of Audi et al. [31]. 

The positron capture rates of a transition from the ith state of the parent to the jth state of the daughter nucleus 
is given by 
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We took the value of D=6295s [32] and the ratio of the axial vector to the vector coupling constant as -1.254 [33]. 
B'^s are the sum of reduced transition probabilities of the Fermi B(F) and GT transitions B(GT). Details of these 
reduced transition probabilities can be found in Ref. [13, 17]. The phase space integral fij is an integral over total 
energy and for positron capture it is given by 
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In above equation w is the total energy of the electron including its rest mass, wi is the total capture threshold energy 
(rest+kinetic) for positron capture. F(-Z,w) are the Fermi functions and were calculated according to the procedure 
adopted by Gove and Martin [34]. G + is the Fermi-Dirac distribution function for positrons. 
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here E = (w — 1) is the kinetic energy of the positrons, Ef is the Fermi energy of the positrons, T is the temperature, 
and k is the Boltzmann constant. 

The number density of electrons associated with protons and nuclei is pY e NA (p is the baryon density, Y e is lepton 
to baryon ratio, and Na is Avogadro's number) 
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here p = (w 2 — I) 1 / 2 is the positron momentum and Eqt. (5) has the units of mol cm~ s . G_ is the Fermi-Dirac 
distribution function for electrons. 
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Eqt. 5 was used for an iterative calculation of Fermi energies for selected values of Y e and T. There is a finite 
probability of occupation of parent excited states in the stellar environment as result of the high temperature prevailing 
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in the interior of massive stars. Weak interactions then also have a finite contribution from these excited states. The 
total positron capture rate per unit time per nucleus is given by 



The summation over all initial and final states was carried out until satisfactory convergence in the rate calculations 
was achieved. Here Pi is the probability of occupation of parent excited states and follows the normal Boltzmann 
distribution. 

The neutrino energy loss rates can occur through four different weak-interaction mediated channels: electron and 
positron emissions, and, continuum electron and positron captures. The neutrino energy loss rates were calculated 
using the same formalism described above except that the phase space integral was replaced by 



For the decay channel Eqt. 8 was used for the calculation of phase space integrals. Upper signs were used for the case 
of electron emissions and lower signs for the case of positron emissions. Regarding the capture channels, Eqt. 9 was 
used for the calculation of phase space integrals keeping upper signs for continuum electron captures and lower signs 
for continuum positron captures. 

The total neutrino energy loss rate per unit time per nucleus is given by 




(7) 
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and by 




(9) 




(10) 



where \\, is the sum of the electron capture and positron decay rates for the transition i — > j. 



On the other hand the total antineutrino energy loss rate per unit time per nucleus is given by 




(11) 



where A-'- is the sum of the positron capture and electron decay rates for the transition i — > j. 



III. UNCERTAINTIES OF THE PN-QRPA MODEL 



The uncertainties involved in stellar rate calculations are considerable. The prevailing extreme physical conditions 
and model parameters invoke uncertainties in the calculation. Lack of experimental data in this scenario deteriorates 
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the situation. The "electron capture direction" can be explored experimentally by (n,p) experiments whereas the "beta 
minus decay direction" can be explored by the (p,n) reactions. There are a handful of other experiments which have 
been used by many theorists to shape the centroid and width of the GT strength. However these experimental data 
are not enough to completely explore the domain of nuclei which are interesting from astrophysical viewpoint. There 
is no experimental data concerning GT strength distribution from parent excited states. In the stellar environment, 
at high temperatures and densities, there is a finite probability of occupation of parent excited states and transitions 
from these excited states are sometimes many orders of magnitude higher than transitions from ground states [19]. 
The weak interaction rates are calculated using 

\ij = -p-fijBij. (12) 

Here i represents the parent excited states and j the daughter's. The first factor is a constant and the second are phase 
space integrals which can be calculated relatively accurate. It is the third factor (the reduced transition probabilities) 
which contains interesting nuclear physics and incorporates uncertainties in the model. Again the reduced transition 
probabilities is a sum of Fermi and GT component (see Eqt. 2). Whereas the calculation of Fermi transition is rather 
straightforward it is precisely the calculation of the excited states and reduced transition probabilities of the GT 
transitions which is the main cause of uncertainty of the underlying model. The pn-QRPA model constructs parent 
and daughter excited states and also calculates GT strength distribution among these states in a microscopic fashion. 
In other words the Brink's hypothesis is not employed in this calculation which increases the reliability of the pn-QRPA 
calculations (Brink's hypothesis states that GT strength distribution on excited states is identical to that from ground 
state, shifted only by the excitation energy of the state) . As mentioned above there are still uncertainties present due to 
the parameters of the model. Roughly, the parameters of the pn-QRPA model can be divided into two different groups: 
(i) "internal" parameters of the model which are by some means adjustable (the pairing gaps and the GT strength 
parameters) and (ii) "external" parameters for which input from other sources like mass formulae (or experimental 
data, if available) is necessary (these include single particle energies and wavefunctions, deformations, Q values and 
neutron/proton separation energies). Whereas "internal" parameters are of minor importance the uncertainty in the 
"external" parameters must be viewed as the limiting factor for the calculation of weak interaction rates of unstable 
nuclide. The values taken for these parameters and their optimal selection procedure were highlighted in the previous 
section. In order to further increase the reliability of the calculated rates experimental data were incorporated into 
the model as also discussed in the previous section. 

We, however, do have a reasonable amount of experimental data on measured half-lives and as such the pn-QRPA 



theory was tested to check the accuracy of the model against the experimentally known half-lives using the same set 
of parameters. The check was performed in both "beta minus decay" and "electron capture" directions. 

In Tables (I) and (II), N denotes the number of experimentally known half- lives shorter than the limit in the second 
column, n is the number (and percentage) of isotopes reproduced under the condition given in the first column, and 
x is the average deviation defined by 
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(Tj°2 is the calculated half-life using the pn-QRPA model and Ty£ is the corresponding measured half-life.) For 
example, the pn-QRPA reproduces 93% of all experimentally known half-lives shorter than 1 minute for /3+/EC 
within a factor of 10 with an average deviation of x = 1.718 and 95% of all known /3~-decaying nuclei with half-lives 
less than a minute are reproduced within a factor of 5 with an average deviation of x = 1.56. It can be seen from the 
tables that the model works better with increasing neutron excess (corresponding to shorter half- lives), that is, with 
increasing distance from stability. This is in agreement with the expectation, since forbidden transitions is neglected 
in the calculation. This is also a promising feature with respect to the prediction of unknown half-lives (specially for 
unstable isotopes), implying that the predictions are made on the basis of a realistic physical model (see also Table I 
and Table J of Ref. [13] for predictive power of pn-QRPA in "electron capture" direction, and Table K of Ref. [13] 
for the predictive power of the model in the "beta minus decay" direction). 

IV. RESULTS AND COMPARISON 

At temperatures pertinent to supernova environment we have a finite probability of occupation of excited states and 
a microscopic calculation of rates from these excited states is desirable. Earlier, Nabi and Sajjad [19] did point to the 
fact that the Brink's hypothesis ( and back resonances for calculation of beta decay rates) is not a good approximation 
to use in stellar rate calculations. Brink's hypothesis states that GT strength distribution on excited states is identical 
to that from ground state, shifted only by the excitation energy of the state whereas the GT back resonances are 
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states reached by the strong GT transitions in the electron capture process built on ground and excited states. The 
luxury of having a huge model space at our disposal allowed us to perform the calculation of electron capture rates 
from 30 excited states of 55 Co. Table III shows the calculated excited states of 55 Co using the pn-QRPA theory. For 
each parent state we calculated the GT strength distribution in a microscopic fashion to around 200 excited states in 
daughter. 

The GT strength distribution (GT±) from ground state and first two excited states of 55 Co was presented earlier 
[21]. We calculated the position of our GT_ (in this direction a neutron is changed into a proton) centroid around 
9.1 MeV for the ground state. For the first two excited states of 55 Co (@ 2.2 MeV and 2.6 MeV), the corresponding 
centroids were placed at 9.2 MeV and 9.6 MeV, respectively. The total GT strength for positron capture from ground 
state of 55 Co was calculated to be around 17.9. Table IV shows the calculated B(GT_) strength values for the ground 
state of 55 Co. The strengths are given up to energy of 10 MeV in daughter nucleus, 55 Ni. Calculated B(GT_) strength 
of magnitude less than 10~ 3 are not included in this table. 

Recently, we presented the extensive calculation of electron capture rates on 55 Co on a fine temperature-density 
scale where we also discussed our results in detail [23]. Here we would like to present a similar calculation for the 
positron capture and the associated (anti)neutrino energy loss rates due to weak-interaction mediated reactions in 
the core of massive stars. 

Fig.l shows four panels depicting our calculated positron capture rates at selected temperature and density domain. 
The upper left panel shows the positron capture rates in low-density region (p[gcm~ 3 ] = 10 ' 5 , 10 1 ' 5 and 10 2 5 ), the 
upper right in medium-low density region (p[gcm~ 3 ] = 10 3 ' 5 , 10 4 5 and 10 5 5 ), the lower left in medium-high density 
region (p{gcm~ 3 ] = 10 6,5 , 10 7 5 and 10 8 5 ) and finally the lower right panel depicts our calculated positron capture 
rates in high density region (p{gcm~ 3 ] = 10 9 5 , 10 10 5 and 10 11 ). The positron capture rates are given in logarithmic 
scales in units of s _1 . T 9 gives the stellar temperature in units of 10 9 K. One should note the order of magnitude 
differences in positron capture rates as the stellar temperature increases. It can be seen from this figure that in the 
low density region the positron capture rates, as a function of stellar temperatures, are more or less superimposed on 
one another. This means that there is no appreciable change in the rates when increasing the density by an order 
of magnitude. We also observe that the positron capture rates are almost the same for the densities in the range 
(10 — 10 6 )g/cm 3 . However as we go from the medium high density region to high density region these rates start to 
'peel off' from one another. Orders of magnitude difference in rates are observed (as a function of density) in high 
density regions. When the densities increase beyond the above stated range a decline in the positron capture rate 
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starts. For a given density the rates increase monotonically with increasing temperatures. 

Fig. 2 and Fig. 3 depict our calculated neutrino and antineutrino energy loss rates due to 55 Co. It is pertinent to 
mention again that the neutrino energy loss rates (depicted in Fig. 2) contain contributions due to electron capture 
and positron decay on 55 Co whereas the antineutrino energy loss rates (Fig. 3) are calculated due to contributions 
from positron capture and electron decay on 55 Co. The energy loss rates are given in logarithmic scales (in units of 
MeV.s -1 ). The figures again consist of four panels depicting the low, medium-low, medium-high and high density 
domains for the core of a massive star. We note the similarity between the positron capture rates (Fig.l) and the 
antineutrino energy loss rates (Fig. 3). The later are slightly enhanced at the corresponding temperature and density. 
On the other hand the neutrino energy loss rates exhibit an entirely different pattern. We note that the neutrino 
energy loss rates are orders of magnitude greater than the corresponding antineutrino energy loss rates. This is an 
expected result (the Q-value of 55 Co for electron capture/positron decay is 3.452 MeV whereas the Q-value of 55 Co 
in the other direction is -8.692 MeV [31]). 

The comparison of electron and positron capture rates on 55 Co with earlier calculations were presented in Ref. [21]. 
Here we present a comparison of the energy loss rates with the earlier calculations. Fig. 4 presents a comparison 
of our calculated neutrino energy loss rates compared with large-scale shell model [12] and FFN [10] calculations. 
The comparison is presented at densities (10 3 , 10 7 , 10 11 )g.cm~ 3 . Compared to large-scale shell model results, our 
calculations lead to a larger energy being carried away by the neutrinos and hence favor cooler cores. We note that 
our corresponding numbers are roughly as big as two orders of magnitude (at presupernova temperature and density 
region). In high temperature and density regions our calculated neutrino energy loss rates are in good comparison 
with those of large-scale shell model. As far as comparison with the pioneering work of FFN is concerned, we note that 
again our rates are enhanced at presupernova temperature-density domain. However at large stellar temperatures and 
densities, FFN neutrino energy loss rates surpass our rates. There are two main reasons for this enhancement of FFN 
rates. Firstly, FFN placed the centroid of the GT strength at too low excitation energies in their compilation of weak 
rates for odd- A and odd-odd nuclei [35]. Secondly, FFN threshold parent excitation energies were not constrained 
and extended well beyond the particle decay channel. At high temperatures contributions from these high excitation 
energies begin to show their cumulative effect. Simulators should take note of our enhanced neutrino energy loss 
rates at the lower temperatures and densities characteristic of the hydrostatic phases of stellar evolution which may 
affect the temperature and the corresponding lepton-to baryon ratio which becomes very important going into stellar 
collapse. 
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The story is different for the comparison of antineutrino energy loss rates (Fig. 5). This time we note that the large- 
scale shell model rates and FFN rates are much more enhanced compared to pn-QRPA rate calculations. Nevertheless 
these are very small numbers and can change by orders of magnitude by a mere change of 0.5 MeV, or less, in parent 
or daughter excitation energies and are more indicative of the uncertainties present in the calculation of the excitation 
energies. 

Fig. 6 shows our summed GT_ strength as a function of excited states in the daughter 55 Ni. We note that almost 
all the strength cumulates up to an energy around 12 MeV in 55 Ni. No appreciable strength is seen in 55 Ni at higher 
excitation energies. 

We finally present our calculated positron capture rates, neutrino and antineutrino energy loss rates on a detailed 
temperature-density grid in Table V. Here Column 1 shows the density in logarithmic scales (in units of gcm~ 3 ), 
Column 2, the stellar temperature in units of 10 9 K , Column 3 gives the calculated positron capture rates in logarithmic 
scales (in units of sec -1 ) at the corresponding temperature and density whereas Column 4 and Column 5 display the 
corresponding neutrino and antineutrino energy loss rates again in logarithmic scales (in units of MeV.sec^ 1 ). All 
logarithms are taken to base 10. Tables of rate calculations presented in earlier compilations (e.g. Ref. [10, 12, 13, 
36, 37]) were not presented on a detail temperature-density grid and at times could lead to erroneous results when 
interpolated. We hope that this table will prove more useful for core-collapse simulators. 

V. CONCLUSIONS 

55 Co is advocated to play a key role amongst the iron- regime nuclide controlling the dynamics of presupernova 
evolution of massive stars. The capture rates and (anti) neutrino energy loss rates on 55 Co are used as nuclear 
physics input parameter for multi-dimensional simulations. Reliable and detailed calculations of these weak-interaction 
mediated rates are desirable for these codes. These parameters may fine tune the final outcome of the neutrino 
transport included multi-dimensional models. 

Here we present, for the first time, an extensive calculation of stellar positron capture rates and the (anti)neutrino 
loss rates for 55 Co on a fine temperature-density scale suitable for simulation codes. According to authors in Ref. 
[11] and Ref. [20], 55 Co is a very important nucleus controlling the events during the pre-collapse phase of iron 
cores of massive stars. The calculated neutrino energy loss rates are around two orders of magnitude enhanced as 
compared to large-scale shell model calculations during the hydrostatic phases of stellar evolution. This may affect 
the temperature, entropy and the lepton-to-baryon ratio which becomes very important going into stellar collapse. 
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We will urge simulators to test run our reported weak interaction rates presented here to check for some interesting 
outcome. We are currently in a phase of extending the present work for other nuclide of astrophysical importance 
and hope to report on the outcome of these calculations in near future. 
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Table (I): The accuracy of the pn-QRPA model compared to experimental data (/3 + /EC, taken from Ref. [27]). TV 
denotes the number of experimentally known half-lives shorter than the limit in the second column, n is the number 
(and percentage) of isotopes reproduced under the condition given in the first column, and x is the average deviation 
defined in the text. 



Conditions 


T^(a) < 


N 


n 


n(%) 


X 


Vx* < 10 


10 6 


894 


706 


79.0 


2.057 




60 


327 


304 


93.0 


1.718 




1 


81 


78 


96.3 


1.848 


Vx, < 2 


10 6 


894 


489 


54.7 


1.363 




60 


327 


245 


74.9 


1.308 




1 


81 


59 


72.8 


1.230 
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Table (II): The accuracy of the pn-QRPA model compared to experimental data , taken from Ref. [26]). N 
denotes the number of experimentally known half-lives shorter than the limit in the second column, n is the number 
(and percentage) of isotopes reproduced under the condition given in the first column, and x is the average deviation 
defined in the text. 



Conditions 


7T/T00 < 


N 


n 


n(%) 


X 


Vx, < 10 


10 6 


654 


472 


72.2 


1.85 ± 1.21 




60 


325 


313 


96.3 


1.67 ± 1.02 




1 


106 


105 


99.1 


1.44 ± 0.40 


Vxi < 5 


10 6 


654 


456 


69.7 


1.68 ± 0.76 




60 


325 


307 


94.5 


1.56 ± 0.66 




1 


106 


105 


99.1 


1.44 ± 0.40 


Vxi < 3 


10 6 


654 


420 


64.2 


1.50 ± 0.46 




60 


325 


295 


90.8 


1.46 ± 0.43 




1 


106 


105 


99.1 


1.44 ± 0.40 


Vxi < 2 


10 6 


654 


369 


56.4 


1.37 ± 0.29 




60 


325 


267 


82.2 


1.36 ± 0.29 




1 


106 


96 


90.6 


1.35 ± 0.27 



Table III: Calculated excited states in parent 55 Co using the pn-QRPA theory in units of MeV. 



0.0 2.17 2.57 2.92 3.08 3.30 3.87 4.10 4.48 4.89 
5.20 5.47 5.68 5.99 6.20 6.50 6.79 7.08 7.33 7.65 
7.88 8.16 8.42 8.67 8.96 9.21 9.48 9.74 9.89 10.00 
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Table IV: Calculated B(GT_) values from ground state in 55 Co. The energy scale refers to excitation energies in 
daughter, 55 Ni. 



Energy(MeV) 


B(GT_) 


Encrgy(McV) B(GT_) 


Energy (Me V) 


B(GT_) 


0.00 


3.06E-01 


6.36 


1.80E-02 


8.36 


5.37E-02 


2.46 


1.40E+00 


6.66 


2.91E-02 


8.62 


2.45E-02 


3.01 


4.48E-01 


6.85 


4.47E-02 


8.73 


3.78E-01 


3.45 


1.04E-01 


7.00 


1.12E-02 


8.90 


3.26E-01 


3.57 


1.80E-03 


7.14 


1.06E-02 


9.11 


2.94E+00 


3.68 


1.79E-01 


7.43 


1.19E-01 


9.26 


1.59E+00 


3.90 


1.99E-01 


7.59 


6.53E-03 


9.42 


1.00E+00 


4.11 


2.25E-02 


7.76 


8.77E-02 


9.60 


2.11E+00 


5.82 


1.93E-02 


7.92 


1.62E-02 


9.87 


4.84E-02 


5.93 


3.89E-03 


8.04 


1.92E-01 


10.03 


7.67E-02 


6.13 


1.58E-02 


8.25 


7.51E-02 
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Table V: Calculated positron capture, neutrino and antineutrino energy loss rates on 55 Co for different selected 
densities and temperatures in stellar matter. log(pYe) has units of g/cm 3 , where p is the baryon density and Y e is the 
ratio of the lepton number to the baryon number. Temperatures (Tg) are measured in 10 9 K. X pc are the positron 
capture rates (sec -1 ). A„ are the total neutrino energy loss rates (MeV.s^ 1 ) due to j3 + decay and electron capture. 
A p are the total antineutrino energy loss rates (MeV.s' 1 ) due to (3~ decay and positron capture. All calculated rates 
are tabulated in logarithmic (to base 10) scale. In the table, -100 means that the rate is smaller than 10~ 100 . 
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logpY c Tg 



logpY c T 9 Ap 



0.5 0.5 

0.5 1 

0.5 1.5 

0.5 2 

0.5 2.5 

0.5 3 

0.5 3.5 

0.5 4 

0.5 4.5 

0.5 5 

0.5 5.5 

0.5 6 

0.5 6.5 

0.5 7 

0.5 7.5 

0.5 8 

0.5 8.5 

0.5 9 

0.5 9.5 

0.5 10 

0.5 15 

0.5 20 

0.5 25 

0.5 30 



0.5 

1 

1.5 

2 

2.5 

3 

3.5 

4 



1 4.5 
1 5 



5.5 

6 

6.5 
7 

7.5 



■95.544 
■48.153 
■32.485 
■24.599 
-19.828 
-16.614 
-14.292 
■12.529 
■11.138 
■10.008 
-9.068 
-8.269 
-7.58 
-6.977 
-6.442 
-5.964 
-5.531 
-5.137 
-4.776 
-4.443 
-2.054 
-0.554 
0.505 

1.3 
■96.024 
■48.155 
-32.485 
■24.599 
■19.827 
-16.614 
■14.292 
■12.529 
■11.138 
-10.008 
-9.067 
-8.269 
-7.58 
-6.976 
-6.442 
-5.963 



-3.375 
-3.376 
-3.366 
-3.305 
-3.143 
-2.884 
-2.585 
-2.287 
-2.006 
-1.745 
-1.502 
-1.276 
-1.063 
-0.862 

-0.67 
-0.486 
-0.309 
-0.138 
0.028 

0.19 
1.569 
2.582 

3.34 
3.931 
-3.375 
-3.376 
-3.366 
-3.305 
-3.142 
-2.884 
-2.584 
-2.286 
-2.005 
-1.744 
-1.501 
-1.275 
-1.062 
-0.861 
-0.669 
-0.485 



-95.515 
-48.102 
-32.41 
-24.498 
-19.699 
-16.457 
-14.107 
-12.315 
-10.896 
-9.738 
-8.77 
-7.944 
-7.229 
-6.601 
-6.042 
-5.539 
-5.084 
-4.668 
-4.286 
-3.932 
-1.368 
0.261 
1.422 
2.301 
-95.995 
-48.103 
-32.409 
-24.497 
-19.698 
-16.456 
-14.106 
-12.314 
-10.894 
-9.736 
-8.768 
-7.943 
-7.228 
-6.599 
-6.04 
-5.538 



9.5 
10 
15 
20 
25 
30 



1.5 0.5 
1.5 1 

1.5 1.5 
1.5 2 

1.5 2.5 
1.5 3 

1.5 3.5 
1.5 4 

1.5 4.5 
1.5 5 

1.5 5.5 
1.5 6 

1.5 6.5 
1.5 7 

1.5 7.5 
1.5 8 

1.5 8.5 
1.5 9 

1.5 9.5 

1.5 10 

1.5 15 

1.5 20 

1.5 25 

1.5 30 



0.5 

1 

1.5 

2 

2.5 
3 

3.5 

4 



-5.531 
-5.137 
-4.776 
-4.442 
-2.053 
-0.553 
0.506 
1.301 
-96.522 
-48.16 
-32.486 
-24.599 
-19.827 
-16.614 
-14.292 
-12.529 
-11.138 
-10.008 
-9.067 
-8.269 
-7.579 
-6.976 
-6.441 
-5.963 
-5.53 
-5.137 
-4.775 
-4.442 
-2.053 
-0.553 
0.506 
1.301 
-97.022 
-48.176 
-32.487 
-24.599 
-19.827 
-16.614 
-14.292 
-12.529 



-0.308 
-0.136 

0.03 
0.191 
1.571 
2.584 
3.342 
3.933 
-3.375 
-3.376 
-3.366 
-3.305 
-3.142 
-2.884 
-2.584 
-2.286 
-2.004 
-1.743 
-1.501 
-1.275 
-1.062 

-0.86 
-0.668 
-0.484 
-0.307 
-0.136 

0.03 
0.191 
1.571 
2.585 
3.343 
3.934 
-3.375 
-3.376 
-3.366 
-3.305 
-3.142 
-2.884 
-2.584 
-2.286 



-5.082 
-4.666 
-4.284 
-3.93 
-1.365 
0.264 
1.425 
2.304 
-96.493 
-48.108 
-32.41 
-24.497 
-19.698 
-16.456 
-14.106 
-12.314 
-10.894 
-9.736 
-8.768 
-7.943 
-7.227 
-6.599 
-6.039 
-5.537 
-5.082 
-4.666 
-4.283 
-3.929 
-1.365 
0.265 
1.426 
2.305 
-96.993 
-48.125 
-32.411 
-24.497 
-19.698 
-16.456 
-14.106 
-12.313 



2.5 
2.5 
2.5 



4.5 

5 

5.5 

6 

6.5 
7 

7.5 

8 

8.5 

9 

9.5 
10 
15 
20 
25 
30 



2.5 0.5 
2.5 1 

2.5 1.5 
2.5 2 

2.5 2.5 
2.5 3 

2.5 3.5 
2.5 4 

2.5 4.5 
2.5 5 

2.5 5.5 



6 

6.5 
7 



2.5 7.5 

2.5 8 

2.5 8.5 

2.5 9 

2.5 9.5 

2.5 10 

2.5 15 

2.5 20 

2.5 25 

2.5 30 



■11.138 
■10.008 
-9.067 
-8.269 
-7.579 
-6.976 
-6.441 
-5.963 

-5.53 
-5.137 
-4.775 
-4.442 
-2.052 
-0.553 

0.506 

1.301 
■97.522 
■48.228 
-32.49 

-24.6 
■19.828 
-16.614 
-14.292 
■12.529 
-11.138 
■10.008 
-9.067 
-8.269 
-7.579 
-6.976 
-6.441 
-5.963 

-5.53 
-5.137 
-4.775 
-4.442 
-2.052 
-0.553 

0.506 

1.302 



-2.004 
-1.743 
-1.501 
-1.274 
-1.062 

-0.86 
-0.668 
-0.484 
-0.307 
-0.136 

0.03 
0.191 
1.571 
2.585 
3.343 
3.934 
-3.375 
-3.375 
-3.366 
-3.305 
-3.142 
-2.883 
-2.584 
-2.285 
-2.004 
-1.743 
-1.501 
-1.274 
-1.062 

-0.86 
-0.668 
-0.484 
-0.307 
-0.136 

0.03 
0.192 
1.572 
2.585 
3.343 
3.935 



-10.894 
-9.736 
-8.768 
-7.942 
-7.227 
-6.598 
-6.039 
-5.537 
-5.081 
-4.665 
-4.283 
-3.929 
-1.364 
0.266 
1.427 
2.306 
-97.493 
-48.177 
-32.414 
-24.498 
-19.698 
-16.456 
-14.106 
-12.314 
-10.894 
-9.736 
-8.768 
-7.942 
-7.227 
-6.598 
-6.039 
-5.537 
-5.081 
-4.665 
-4.283 
-3.929 
-1.364 
0.266 
1.427 
2.306 
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logpY e 


Tg 




K 


Ac 


logpY e 


Tg 




K 


Ac 


logpY e 


Tg 


Ape 


K 


Ac 


3 


0.5 


-98.023 


-3.374 


-97.994 


4 


0.5 


-99.036 


-3.364 


-99.006 


5 


0.5 


-100 


-3.281 


-100 


3 


1 


-48.382 


-3.375 


-48.331 


4 


1 


-49.204 


-3.366 


-49.153 


5 


1 


-50.236 


-3.285 


-50.185 


3 


1.5 


-32.501 


-3.366 


-32.425 


4 


1.5 


-32.643 


-3.36 


-32.567 


5 


1.5 


-33.377 


-3.276 


-33.301 


3 


2 


-24.603 


-3.305 


-24.501 


4 


2 


-24.634 


-3.299 


-24.532 


5 


2 


-24.923 


-3.229 


-24.821 


3 


2.5 


-19.828 


-3.142 


-19.699 


4 


2.5 


-19.84 


-3.137 


-19.71 


5 


2.5 


-19.953 


-3.084 


-19.823 


3 


3 


-16.614 


-2.883 


-16.457 


4 


3 


-16.62 


-2.88 


-16.462 


5 


3 


-16.67 


3 -2.843 


-16.515 


3 


3.5 


-14.292 


-2.584 


-14.106 


4 


3.5 


-14.295 


-2.581 


-14.109 


5 


3.5 


-14.325 


-2.558 


-14.138 


3 


4 


-12.529 


-2.285 


-12.314 


4 


4 


-12.53 


-2.284 


-12.315 


5 


4 


-12.549 


-2.269 


-12.333 


3 


4.5 


-11.138 


-2.004 


-10.894 


4 


4.5 


-11.139 


-2.003 


-10.895 


5 


4.5 


-11.151 


-1.993 


-10.907 


3 


5 


-10.008 


-1.743 


-9.736 


4 


5 


-10.008 


-1.742 


-9.737 


5 


5 


-10.017 


-1.735 


-9.745 


3 


5.5 


-9.067 


-1.501 


-8.768 


4 


5.5 


-9.068 


-1.5 


-8.768 


5 


5.5 


-9.074 


-1.495 


-8.774 


3 


6 


-8.269 


-1.274 


-7.942 


4 


6 


-8.269 


-1.274 


-7.943 


5 


6 


-8.274 


-1.27 


-7.947 


3 


6.5 


-7.579 


-1.061 


-7.227 


4 


6.5 


-7.58 


-1.061 


-7.227 


5 


6.5 


-7.583 


-1.058 


-7.231 


3 


7 


-6.976 


-0.86 


-6.598 


4 


7 


-6.976 


-0.86 


-6.598 


5 


7 


-6.979 


-0.857 


-6.601 


3 


7.5 


-6.441 


-0.668 


-6.039 


4 


7.5 


-6.442 


-0.668 


-6.039 


5 


7.5 


-6.444 


-0.666 


-6.041 


3 


8 


-5.963 


-0.484 


-5.537 


4 


8 


-5.963 


-0.484 


-5.537 


5 


8 


-5.965 


-0.482 


-5.539 


3 


8.5 


-5.53 


-0.307 


-5.081 


4 


8.5 


-5.53 


-0.307 


-5.081 


5 


8.5 


-5.532 


-0.305 


-5.083 


3 


9 


-5.136 


-0.136 


-4.665 


4 


9 


-5.137 


-0.135 


-4.665 


5 


9 


-5.138 


-0.134 


-4.666 


3 


9.5 


-4.775 


0.03 


-4.283 


4 


9.5 


-4.775 


0.031 


-4.283 


5 


9.5 


-4.776 


0.032 


-4.284 


3 


10 


-4.442 


0.192 


-3.929 


4 


10 


-4.442 


0.192 


-3.929 


5 


10 


-4.443 


0.193 


-3.929 


3 


15 


-2.052 


1.572 


-1.364 


4 


15 


-2.052 


1.572 


-1.364 


5 


15 


-2.053 


1.572 


-1.364 


3 


20 


-0.552 


2.585 


0.266 


4 


20 


-0.552 


2.585 


0.266 


5 


20 


-0.553 


2.586 


0.266 


3 


25 


0.506 


3.343 


1.427 


4 


25 


0.507 


3.343 


1.427 


5 


25 


0.506 


3.343 


1.427 


3 


30 


1.302 


3.935 


2.306 


4 


30 


1.302 


3.935 


2.306 


5 


30 


1.302 


3.935 


2.306 


3.5 


0.5 


-98.526 


-3.372 


-98.497 


4.5 


0.5 


-99.565 


-3.342 


-99.536 


5.5 


0.5 


-100 


-3.144 


-100 


3.5 


1 


-48.73 


-3.373 


-48.679 


4.5 


1 


-49.71 


-3.345 


-49.659 


5.5 


1 


-50.815 


-3.14 


-50.764 


3.5 


1.5 


-32.536 


-3.365 


-32.46 


4.5 


1.5 


-32.921 


-3.341 


-32.845 


5.5 


1.5 


-33.907 


-3.111 


-33.831 


3.5 


2 


-24.61 


-3.303 


-24.508 


4.5 


2 


-24.709 
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FIG. 1: Positron capture rates on 55 Co, as a function of stellar temperatures, for different selected densities . Densities are 
in units of gcm~ 3 . Temperatures are measured in 10 9 K and logioA pc represents the log of positron capture rates in units of 
sec -1 . 
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FIG. 3: same as Fig. 2 but for antineutrino energy loss rates. 
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FIG. 4: Comparison of neutrino energy loss rates with those of large-scale shell model [12] and FFN [10] calculations as a 
function of stellar temperatures and densities. 




FIG. 5: Same as Fig.4 but for antineutrino energy loss rates. 




FIG. 6: Cumulative sum of the B(GT_) values. The energy scale refers to excitation energies in daughter 55 Ni. 
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